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The nature and origin of turbulence and magnetic fields in the intergalactic 
space are important problems that are yet to be understood. We propose a 
scenario in which turbulent flow motions are induced via the cascade of the 
vorticity generated at cosmological shocks during the formation of the large 
scale structure. The turbulence in turn amplifies weak seed magnetic fields of 
any origin. Supercomputer simulations show that the turbulence is subsonic 
inside clusters/groups of galaxies, whereas it is transonic or mildly supersonic 
in filaments. Based on a turbulence dynamo model, we then estimate that the 
average magnetic field strength would be a few fiG inside clusters/groups, 
approximately 0.1 /iG around clusters/groups, and approximately 10 nG in 
filaments. Our model presents a physical mechanism that transfers the grav- 
itation energy to the turbulence and magnetic field energies in the large scale 
structure of the universe. 
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There is growing evidence that the intergalactic medium (IGM) is permeated with mag- 
netic fields and is in a state of turbulence, similar to the interstellar medium within galax- 
ies. Magnetic fields in the intracluster medium (ICM) have been measured using a variety of 
techniques, including observations of diffuse synchrotron emission from radio halos, inverse- 
Compton scattered cosmic background radiation in extreme ultraviolet and hard x-ray radiation, 
and Faraday rotation measure (RM). The inferred strength of the magnetic fields is on the order 
of 1 nG CDI2I2]). In the IGM outside of clusters, an upper limit of ~ 0.1 fiG has been placed on 
the magnetic field strength of filaments, based on the observed limit of the RMs of background 
quasars d^JEl)- 

So far, signatures of turbulence have been observed only in the ICM. The analysis of the 
gas pressure maps of the Coma cluster revealed that pressure fluctuations are consistent with 
Kolmogoroff turbulence, and turbulence is likely to be subsonic with e tuJ . h £ 0.1e t h» where 
£ tur b and e t h are the turbulence and thermal energy densities, respectively ©. These results 
agree with predictions of numerical simulations of large scale structure (LSS) formation (]7]IU>. 
Turbulence in the ICM also has been studied in RM maps of a few clusters (1911701). 

It has been suggested that cosmological shocks with Mach numbers up to ~ 10 4 and speeds 
up to a few thousand km s _1 exist in the IGM ([77] [72] [721). Such shocks result from the super- 
sonic flow motions that are induced by the hierarchical formation of LSS in the universe. They 
are collisionless shocks, which form in a tenuous plasma via collective electromagnetic inter- 
actions between particles and electromagnetic fields ((7?1) . The gravitational energy released 
during the structure formation is transferred by these shocks to the IGM plasma in several 
different forms: in addition to the gas entropy, cosmic rays are produced via diffusive shock 
acceleration 1/5117(51) . magnetic fields are generated via the Biermann battery mechanism (1771 [71) 
and Weibel instability jl8\\19\i . and vorticity is generated at curved shocks ([201127]) . 

In astrophysical plasmas in which charged particles are coupled to magnetic fields, turbulent 
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flow motions and magnetic fields are closely related. We suggest that the turbulence in the IGM 
is induced by the cascade of the vorticity generated at cosmological shocks. The turbulence 
then amplifies the intergalactic magnetic fields (IGMFs) through the stretching of field lines, 
a process known as the turbulence dynamo. This scenario provides a theoretically motivated 
model for the evolution of the IGMFs in LSS, independent of the origin of seed fields. 

There are other sources that can also provide turbulence and magnetic fields to the IGM. For 
instance, galactic winds can drag out the galactic magnetic fields on the order of 1 /xG strength 
into the surrounding IGM (|22j) . The magnetic fields in the lobes of the jets from galactic black 
holes can also contaminate the IGM (23). Mergers of smaller objects are expected to produce 
turbulent motions in the ICM, which in turn amplify the existing magnetic fields (|2?1) . Those 
processes, although possibly important, are not topics of this study. 

We first calculated the vorticity, uj = V x v (curl of flow velocity), in the IGM, from a 
numerical simulation using particle-mesh/Eulerian hydrodynamic code d25b for the formation 
of LSS in a cold dark matter dominated universe with a cosmological constant [supporting 
online material (SOM) text SI]. As shown in Fig. 1, numerous shocks exist in the LSS that are 
bounded by accretion shocks (11). The distribution of vorticity closely matches that of shocks, 
suggesting that a substantial portion of the vorticity, if not all, must have been generated at the 
shocks. 

There is a clear trend that the vorticity is larger in hotter (Fig. 2) and denser (fig. SI) regions. 
As shown in the top right panel of Fig. 2, at the present epoch, tu rms t a ge ~ 10 to 30 (uj vms , 
the root mean square of the vorticity; t age , the present age of the universe) in clusters/groups 
(temperature T > 10 7 K) and filaments (10 5 < T < 10 7 K), whereas it is on the order of unity 
in sheetlike structures (10 4 < T < 10 5 K) and even smaller in voids (T < 10 4 K) (see SOM text 
S2 for the temperature phases of the IGM). It increases a little with time and asymptotes after 
redshift z £ L Because the local eddy turnover time, t cddy , can be defined with the vorticity as 
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teddy — 1/uj, u>t age (z) represents the number of eddy turnovers in the age of the universe at a 
given z. Roughly, if tut^ is greater than a few, we expect there has been enough time for the 
vorticity to cascade down to smaller scales and for turbulence to develop in the IGM. So it is 
likely that turbulence is well developed in clusters/groups and filaments, but the flow is mostly 
non-turbulent in sheets and voids. 

In our simulation the vorticity was generated either directly at curved cosmological shocks 
or by the baroclinity of flows. The baroclinity resulted from the entropy variation induced at 
shocks. Therefore, the baroclinic vorticity generation also can be attributed to the presence of 
cosmological shocks. Our estimates of vorticity generation by the two processes (SOM text S3) 
are shown with open symbols in the top right panel of Fig. 2. They agree reasonably well with 
the vorticity present in the simulation, although the estimates are intended to be rough. The plot 
indicates that the contributions from the two processes are comparable. 

To estimate the energy associated with turbulence, the curl component of flow motions, ?; cur i, 
which satisfies the relation V x v CUI \ = V x v, is extracted from the velocity field (SOM text 
S4). As vorticity cascades to develop into turbulence, the energy (l/2)pv^ ml (p, gas density) is 
transferred to turbulent motions, so we regard it as the turbulence energy, £turb- As shown in Fig. 
3, £turb < £th m clusters/groups. In particular, the mass-averaged value is (£turb/£th) mass = 0.1 
to 0.3 for T > 10 7 K, which is in good agreement with the observationally inferred value in 
cluster cores ©. The turbulence Mach number M turb = v tuvh /c s = vT8 (^turb/^th) 1 ^ 2 , where 
c s is the sound speed. Therefore, overall turbulence is subsonic in clusters/groups, whereas it is 
transonic or mildly supersonic in filaments. 

The general consensus regarding the origin of the IGMFs is that no mechanism can produce 
strong coherent magnetic fields in the IGM before the formation of LSS and galaxies (l2oT) . 
However, it is reasonable to assume that weak seed fields were created in the early universe 
(SOM text S5). The seed fields can be amplified by the intergalactic turbulence discussed above. 
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In principle, if we were to perform magnetohydrodynamic (MHD) simulations of structure 
formation, the amplification of the IGMFs could be followed. In practice, however, the currently 
available computational resources do not allow a numerical resolution high enough to reproduce 
the full development of MHD turbulence in LSS ©. 

In order to follow the growth of the IGMFs by the dynamo action of turbulence, we turned 
to a separate simulation in a controlled box. Starting with a very weak regular field, a three- 
dimensional incompressible simulation of driven MHD turbulence was performed (SOM text 
S6). In the simulation, the evolution of magnetic fields goes through three stages: (i) the initial 
exponential growth stage, when the back-reaction of magnetic fields is negligible; (ii) the linear 
growth stage, when the back-reaction starts to operate; and (iii) the final saturation stage (|27l) . 
Adopting the simulation result, we model the growth and saturation of magnetic energy as 

[ 0.04 x exp [(t/tcdd y - 4)/0.36] for t/t cddy < 4 
0(t/tcdd y ) = — = I (0.36/41) x (t/t eddy - 4) + 0.04 for 4 < t/t eddy < 45 (1) 
£turb [0.4 for t/t cddy > 45 

(fig. S2). Assuming that the fraction of turbulence energy governed by Eq. 1, 0, is converted 
into the magnetic energy, we estimate the strength of the IGMFs as B = [87T£turb • (K^agc)] 1 ^ 2 - 
Here the values of u and £ tU rb are calculated locally from the structure formation simulation. 

The resulting IGMFs follows the cosmic web of matter distribution as shown in Fig. 4 
(and in fig. S3). On average the IGMFs are stronger in hotter (Fig. 2) and denser (fig. SI) 
regions in our model. The strength of the IGMFs is B £ 1 jiG inside clusters/groups (the mass- 
averaged value for T > 10 7 K), ~ 0.1 /xG around clusters/groups (the volume- averaged value 
for T > 10 7 K), and ~ 10 nG in filaments at present (bottom right panel of Fig. 2) (see SOM 
text S7 for the numerical convergence of the estimation). These values agree with the observed 
field strengths discussed earlier. They also agree with the previous study ©, in which the 
magnetic field strength in clusters was estimated to be a few /iG, based on a kinetic theory. The 
IGMFs should be much weaker in sheetlike structures and voids. But as noted above, turbulence 
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is not fully developed in such low density regions, so our model is not adequate to predict the 
field strength there. For each temperature phase, the IGMFs are stronger in the past, because 
the gas density is higher. However, the IGMFs averaged over the entire computational volume 
are weaker in the past because the fraction of strong field regions is smaller. 

While being amplified, magnetic fields become coherent through the inverse cascade (l27l) . 
The coherence scale of magnetic fields in fully developed turbulence is expected to be several 
times smaller than the driving scale, that is, the scale of dominant eddies (SOM text S8). In the 
IGM outside of clusters, the curvature radius of typical cosmological shocks is approximately a 
couple of Mpc (177]) (fig. S4), which should represent a characteristic scale of dominant eddies. 
The coherence length of the IGMFs there is then expected to be several 100 kpc. On the other 
hand, the scale height of the ICM is several 100 kpc. If it corresponds to the scale of the 
dominant eddies, the coherence length in the ICM is expected to be ~ 100 kpc or so. 

Our model can predict the RMs owing to the IGMFs, which may be tested in future observa- 
tions with Low Frequency Array and Square Kilometer Array (|28|) . Also, our model IGMFs can 
be employed in the study of the propagation of ultra-high-energy cosmic rays, which is crucial 
to search for astrophysical accelerators of such high energy particles (|29|) . 
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Fig. 1. Two-dimensional images showing gas density p in a logarithmic scale (left), locations of 
shocks with color-coded shock speed f s hock (middle), and magnitude of vorticity cut agc (right), 
around a cluster complex of (25 h^Mpc) 2 area at present (z = 0). Here, h is the Hubble 
constant in units of 100 km s _1 Mpc _1 . The complex includes a cluster of x-ray emission- 
weighted temperature T x m 3.3 keV. Color codes for each panel are (left) pj (p) from 10 _1 
(green) to 10 4 (red); (middle) f s h oc k from 15 (green) to 1,800 km s -1 (red); and (right) a>t ag e 
from 0.5 (green) to 100 (red). 
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Fig. 2. (Left) Volume fraction with given temperature and vorticity magnitude (top left) and 
temperature and magnetic field strength (bottom left) at present. (Right) Time evolution of the 
root mean square of the vorticity (top right) and the volume-averaged magnetic field strength 
(bottom right) for four temperature phases of the IGM and for all the gas as a function of redshift 
z. Magenta symbols in the top right panel are our estimates of the vorticity generated directly 
at curved shocks (open circles) and by the baroclinity of flows (open squares). Magenta open 
circles in the bottom right panel show the mass-averaged magnetic field strength for T > 10 7 
K. 
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Fig. 3. Ratio of turbulence to thermal energies as a function of temperature at present. The 
values shown are volume- averaged and mass-averaged over temperature bins. 
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Fig. 4. Volume-rendering image showing the logarithmically scaled magnetic field strength at 
z = in the whole computational box of (100 /i _1 Mpc) 3 volume. Color codes the magnetic 
field strength from 0. 1 nG (yellow) to 10 fiG (magenta). The colors were chosen so that clusters 
and groups show as magenta and blue and filaments as green. 
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Supporting Online Material 
Materials and Methods 

51. Simulation of large scale structure formation 

The simulation used to estimate the vorticity in the intergalactic medium (IGM) is for a ACDM 
universe with the following values of cosmological parameters: VL B m = 0.043, Vt DM = 0.227, 
Vt K = 0.73, h = i7 /(100km/s/Mpc) = 0.7,ando- 8 = 0.8. It was performed using a PM/Eulerian 
hydrodynamic code for the formation of large scale structure (LSS) (SI), in a cubic region of 
comoving volume (100 /i _1 Mpc) 3 with 1024 3 grid zones for gas and gravity and 512 3 particles 
for dark matter, allowing a spatial resolution of Al = 97.7 kpc. The simulation is adiabatic 
in the sense that it does not include radiative cooling, galaxy/star formation, feedbacks from 
galaxies/stars, and reionization of the IGM. Although these processes play important roles in 
determining the physical conditions mainly in cluster cores and voids, it was shown that they do 
not affect significantly the global properties and dissipations of cosmological shocks (52, S3). 

52. Phases of the IGM 

The intergalactic gas is heated mainly by cosmological shocks (S4, S5). It was suggested that 
the IGM can be classified roughly into four phases according to gas temperature: the hot gas 
with T > 10 7 K mostly inside and around clusters/groups, the warm-hot intergalactic medium 
(WHIM) with T = 10 5 — 10 7 K found mostly in filaments, the low temperature WHIM with 
T = 10 4 — 10 5 K distributed mostly as sheetlike structures, and the diffuse gas with T < 10 4 K 
residing mostly in voids (S5, S6). In particular, the intracluster medium (ICM) refers to the 
hottest phase of the IGM inside clusters that is observable in x-ray. 

53. Vorticity generation 

Vorticity can be generated directly at curved shocks and by the baroclinity of flows. For uniform 
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upstream flow, the vorticity produced behind curved shock surface is 

£ cs = ^ p2 ~ Px)2 KXJ X x n, (SI) 
P1P2 

— * 

where p\ and p 2 are the upstream and downstream gas densities, respectively, U\ is the upstream 
flow velocity in the shock rest frame, K is the curvature tensor of the shock surface, and n is 
the unit vector normal to the surface (57). If isopycnic surfaces (constant density surfaces) do 
not coincide with isobaric surfaces, vorticity is generated with the rate given by 

c^bc = \v P x Vp. (S2) 
P 

Rough estimates of the vorticity generations by the two processes in the simulation of LSS 
formation are made as follows. Considering the generation is a random walk process, 

(wAt) m (^^) V2 (53) 

is calculated. Here At is the duration of coherent generation, and so t agc (z)/At is the number 
of random walks during the age of the universe at redshift z. For the direct generation at 
cosmological shocks, 

{p2 - pi) 2 U nl U n2 g 
P1P2 R Al 

is used. Here R is the curvature radius of the shocks, U n i and U n2 are the flow speeds normal 
to the shocks in upstream and downstream, Al is the grid size, and g takes account of all other 
factors such as geometric ones. For the baroclinic vorticity generation, uj hc from Eq. (52) is 
used. In the estimations, the distribution and property of shocks and flows from the simulation 
are used. Also At is taken to be 10 8 yrs, which is a typical dynamical time, e.g., the sound 
crossing time in cluster cores. The comoving value of (R) ~ 1 /i _1 Mpc/(l + z), and (g) ~ 1 
are taken. Note that 1 h^Mpc is the typical scale of nonlinear structures at present. The factor 
1/(1 + z) is included because the average curvature radius of cosmological shocks is smaller in 
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the past even in the comoving frame (S8). In Fig. 2 the quantities normalized with l/t age (z) are 
shown. 

54. Decomposition of flow velocity 

The flow velocity can be decomposed into 

v = v div + V CUT \ + tf uni f , (S5) 

where the divergence and curl components are defined as V ■ v div = V • v and V x v cm i eVxw, 
respectively. That is, v div is associated to compressional motions, while v cm i to incompressible 
shear motions. Here v unii is the component uniform across the computational box, whose mag- 
nitude is much smaller than the other two components. The decomposition is calculated exactly 
in Fourier space. We note with the above decomposition, locally 

Vdiv • ^curl ^0 SO -V 2 ^ - (vj iv + V 2 UTl ) . (S6) 

However, globally 

/ ^div • v CUTl d 3 x = so f l-v 2 d 3 x = I - (vl iv + v 2 url ) d 3 x. (S7) 

J box J box Z J box Z v ' 

55. Seed magnetic fields 

A number of mechanisms that can create seeds field in the early universe have been suggested. 
Besides various inflationary and string theory mechanisms, the most popular astrophysical 
mechanism is the Biermann battery (59). It was discussed in the context of cosmological shocks 
(S10) and ionization fronts (Sll), and shown to build up weak fields of strength up to 1CT 20 G 
by z ~ a few. 

At cosmological shocks, in addition, Weibel instability can operate and produce magnetic 
fields up to the level of eb ~ lCT 3 £ s h (S12), and streaming cosmic rays accelerated by the 
shocks can amplify weak upstream magnetic fields up to the level of e B ~ (l/2)(L r i/c)e C Rs 
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via non-resonant growing mode (SI 3). Here £ s h is the energy density of upstream flow, and Eb 
and £cr s are the energy densities of downstream magnetic fields and CRs. With U\/c ~ 1CT 3 
for cosmological shocks (S2, S3), these processes can potentially produce stronger seed fields, 
although the coherence length of the resulting fields is expected to very small and microscopic. 

S6. Simulation of MHD turbulence 

First we note there are at least two distinct types of magnetohydrodynamic (MHD) turbulence: 
one with strong regular fields and the other with weak/zero regular fields. The former has 
been described successfully with the nonlinear interactions of Alfven waves (S14). The latter is 
mainly hydrodynamic in large scales, but is modified by dynamically important magnetic fields 
in small scales (S15, S16). It is the latter that is relevant to this work. 

Incompressible, driven turbulence was simulated initially with a very weak regular field. 
A pseudo-spectral code (S15) was used, employing hyper-viscosity and hyper-resistivity with 
the Prandtl number of unity. The advantages of this numerical approach include virtually zero 
intrinsic numerical viscosity and resistivity, and the maximized inertial range. The simulation 
was performed with a numerical resolution of 256 3 collocation points. The turbulence was 
driven at the scale of L driv m g ~ (l/2)L box , where L hox is the computational box size. The 
driving strength was set so that the total turbulence energy becomes £ tU rb = £ km + £b ~ 1 at 
saturation. Initially e B = 1CT 6 , but the evolution is not sensitive to the initial field strength 
as long as it is sufficiently weak (S16). The left panel of Fig. S2 shows the time evolution 
of kinetic and magnetic energies. Here, the eddy turnover time is defined with the vorticity at 
driving scale at saturation, t cddy = l/c<j driving . 

We made the incompressible simulation, because compressible simulations need much higher 
resolution to achieve the same growth rate and magnetic field strength at saturation. In addi- 
tion, due to large numerical dissipation, controlling viscosity and resistivity is not trivial in 
compressible simulations. Nevertheless, in order to show the general behavior of turbulence 
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dynamo seen in our incompressible simulation remains the same in compressible regime, we 
performed compressible simulations of driven MHD turbulence, using a code based on the 
third-order essentially non-oscillatory (ENO) upwind scheme {SIT). This is one of the numeri- 
cal schemes for MHD with least numerical dissipation. Different numerical resolutions of 8 3 to 
216 3 grid zones were used. But otherwise, the initial conditions, the driving of turbulence, and 
hyper-viscosity and hyper-resistivity with the Prandtl number of unity were the same as in the 
incompressible simulation. We only considered the case with the Mach number M s ~ 1, which 
is most relevant to this work. We note that turbulence dynamo can be suppressed in highly 
supersonic flows (S18). 

Fig. S5 compares the growth rate (left panel) and the magnetic energy at saturation (right 
panel) from compressible simulations with those from incompressible simulation. Two points 
are clear from the figure: (1) the growth rate is slower in compressible simulations, but the 
pattern follows that in incompressible turbulence, and (2) compressible simulation would need 
more than 1000 3 grid zones to achieve the saturated magnetic field strength of incompressible 
simulation. It also demonstrates why it is impossible to reproduce the full development of 
MHD turbulence in the simulations of LSS formation with the currently available computational 
resources. 

S7. Convergence of our results 

In order to test the numerical convergence of our estimation for the strength of the intergalactic 
magnetic fields (IGMFs), we repeated the same analysis for the simulations of LSS formation 
with lower numerical resolutions of 64 3 to 512 3 grid zones. Except the resolution, they are the 
same simulations as that in §S1, in the sense that they have the same box size, the same realiza- 
tion of initial conditions, and the same physics included. Fig. S6 shows the volume- averaged 
magnetic field strength with different resolutions. The averaged strength around clusters/groups 
with T > 10 7 K is well converged at the resolution of 1024 3 . On the other hand the averaged 

17 



strength in filaments with 10 5 < T < 10 7 K would be underestimated. The converged value 
may be ~ 2 — 3 times larger than our estimate with the resolution of 1024 3 . 

S8. Coherence length 

The consideration of the coherence length of turbulent magnetic fields raises the issues of the 
scale of energy equipartition and the magnetic field structure at the scale. The right panel of Fig. 
S2 shows the power spectra for flow velocity and magnetic fields at a time of saturation from the 
simulation in §S6. The energy equipartition occurs at ~ (l/2.5)L driving , which is close to the 
scale of the magnetic energy peak, ~ (l/3)L driv ing- The kinetic energy peak occurs at the the 
driving scale, L dr i ving . Above the equipartition scale the flow structure is nearly isotropic, while 
below it both the flow and magnetic field structures are anisotropic and the eddies are stretched 
along the local magnetic field lines (S15). On the other hand, the magnetic field has most power 
at ~ (l/6)L driving . So we may argue that in the fully developed stage of the MHD turbulence 
considered here, the coherence length of magnetic fields would be several times smaller than 
the driving scale, the scale of dominant eddies. 

Additional Figures Used in Manuscript 

Fig. SI shows the correlation between the magnitude of vorticity and the gas density (left panel) 
and the strength of magnetic fields and the gas density (right panel) at present. 

Fig. S3 shows the spatial distribution of our model IGMFs around a cluster complex and along 
a filament. In the left image, in addition to the main cluster of T x rs 3.3 keV, a group has B ~ a 
few fiG. They are surrounded by a broad region of B ~ 10 nG. In the right image, groups with 
B ~ jiG are distributed along the filament with B ~ 10 nG. 

Fig. S4 shows the power spectra for the flow velocity and its curl and divergence components 
(§S4) at present from the simulation of LSS formation in §S1. At long wavelengths, the ampli- 
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tude of perturbations are small, so that the linear theory applies. That is, P cm i(k) — > as k — > 0, 
while Pdiv(k) follows the analytic theory expectation, Pdiv(^) ~ k^ 1 . For wavelengths smaller 
than a few /i _1 Mpc, nonlinearities dominate, and we see that P CU ri(^) ~ -Pdiv(^)- P C ui\(k) 
peaks at ~ 5 h^Mpc, and for values of k somewhat larger than the peak wavenumber, the 
spectrum follows a power law of k~ 5 ^ 3 , the Kolmogorov spectrum. P cur \(k) has most power at 
~ 2 h^Mpc, that indicates the typical scale of nonlinear structures in the simulation. 
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Fig. SI. Volume fraction with given gas density and vorticity magnitude (left panel) and gas 
density and magnetic field strength (right panel) at present. 
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Fig. S2. Left panel: Time evolution of kinetic and magnetic energies in a three-dimensional, 
incompressible simulation of driven MHD turbulence with very weak initial magnetic field. 
The green lines show our fitting for the growth and saturation of magnetic energy. Right panel: 
Power spectra for flow velocity, P v , and magnetic fields, P B , at a time of saturation. Two 
straight lines of slopes —5/3 and —1 are also drawn for comparison. 



21 




Fig. S3. Volume rendering images showing logarithmically scaled magnetic field strength at 
z = in a volume of (25 /i _1 Mpc) 3 around the same cluster complex as shown in Fig. 1 (left 
panel) and in a volume of 25 x 15.6 x 6.25 (/i _1 Mpc) 3 which includes a number of groups 
along a filament (right panel). As in Fig. 4, color codes the magnetic field strength from 0.1 nG 
(yellow) to 10 fj,G (magenta). Clusters and groups are shown with magenta and blue while 
filaments with green. 



22 




log(k/27r) (l/lOOh^Mpc) 

Fig. S4. Power spectra, / P v dk = ((l/2)t> 2 ), for flow velocity and its curl and divergence 
components at present in the simulation of LSS formation. Two straight lines of slopes —5/3 
and —1 are also drawn for comparison. 
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Fig. S5. Left panel: Time evolution of magnetic energy in three-dimensional, compressible 
simulations of driven MHD turbulence with very weak initial magnetic field. For comparison, 
the evolution in incompressible simulation shown in Fig. S3 is also plotted. Right panel: Mag- 
netic energy at saturation in compressible simulations with different resolutions of A^ u . The 
value in incompressible simulation is marked with solid line. 
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Fig. S6. Volume- averaged magnetic field strength at present for four temperature phases of the 
IGM and for all the gas from structure formation simulations with different resolutions of iV c 3 e ii- 
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